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Abstract 

We have studied the position dependence of neutrino energy on the Kusenko- 
Segre mechanism as an explanation of the proper motion of pulsars. The 
mechanism is also examined in three-generation mixing of neutrinos and in a 
non-adiabatic case. The position dependence of neutrino energy requires the 
higher value of magnetic field such as B ~ 3 x 10 15 Gauss in order to explain 
the observed proper motion of pulsars. It is shown that possible non-adiabatic 
processes decrease the neutrino momentum asymmetry, whereas an excess of 
electron neutrino flux over other flavor neutrino fluxes increases the neutrino 
momentum asymmetry. It is also shown that a general treatment with all 
three neutrinos does not modify the result of the two generation treatment if 
the standard neutrino mass hierarchy is assumed. 
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In a recent paper, Kusenko and Segre [|IJ have proposed a new explanation for the peculiar 
velocities (or kick velocities) of pulsars by considering magnetically distorted neutrino oscil- 
lations. Many mechanisms for the phenomena have been suggested by some asymmetry 
during the supernova collapsing and explosion. However, most of them have difficulties in 
explaining the observed average pulsar velocity 450 ± 90 km/s [Q. On the other hand, 
the mechanism suggested in Ref. Q, hereafter called the KS mechanism, is successful |7| 
for a reasonable value of the magnetic field in pulsars, i.e., B ~ io 14 ~ 15 Q auss The KS 
mechanism assumes that neutrino oscillations, v e «-> u T , take place between the electron- 
and tau-neutrinospheres, where the matter density is p ~ io n ~ 12 g/cm 3 . Since neutrinos 
emitted during the cooling of a protoneutron star have roughly 100 times the momentum 
of the proper motion of the pulsar, the neutrino momentum asymmetry of Ak/k ~ 1 % 
can explain the observed motion of pulsars. Reference shows that the magnetic field of 
3 x 10 14 Gauss can produce the desirable momentum asymmetry. In his comment, Qian || 
pointed out that the approximate relation used in Ref. |I| is not a valid assumption and that 
somewhat higher magnetic field is required in order to explain the observed kick velocity. 
The differences between Ref. and Ref. || are in the flux factor 1/3 and the estimate of 
the matter density (p = 10 11 and 10 12 g/cm 3 ) besides the different expressions for the ratio 
of scale heights h Ne /h T , where h Ne = |d In A^ e /dr| _1 and h T = |dlnT/dr| _1 . 

Therefore, it is worthwhile to re-examine the KS mechanism in a more general case for 
supernovae and neutrinos. In this paper, first we will consider the position dependence of 
neutrino energy inside the protoneutron star, which was ignored in both Ref. ||] and 0. Our 
expression for the neutrino momentum asymmetry will be shown to reproduce the results 
in Ref. || and Ref. U, respectively, under different approximations. We will also study 
the treatment of the KS mechanism in a three-neutrino generation scheme and see how the 
results are modified for possible non-adiabatic processes. The possibility of different fluxes 
for each neutrino flavor and its effect on the pulsar velocities are also discussed. 

First, let us briefly review the KS mechanism in which only v e <-> v T conversions are 
considered. The difference in the mean free paths of v e and v T renders the v e neutrinosphere, 
S e , to be located at larger radius than the v T neutrinosphere, S T . If a resonant conversion 
between v e and v T occurs outside S T but inside S e , the v e converted from v T is absorbed 
or thermalized by the background medium whereas the v T converted from v e escapes freely. 
In effect, the z/ r 's are expected to propagate out not from the S T , but from the resonance 
surface which is distorted in the presence of magnetic field. The resonance condition for two 
neutrino oscillations in B field is || 

A ° cos20 = V2G F N e (r) + ^[^^\ k-B (1) 




2k(r 

= V2G F N t 
= V2G F N e [l + e B ], 

where Aq = m{v^) 2 — m{v\) 2 in the case of v e <-> v T conversion, N e = (n e - — n e +), 9 
is the neutrino mixing angle in vacuum, Gp the Fermi constant, B the magnetic field of 
pulsars, k ~ \k\ the energy of relativistic neutrinos, 4> the angle between B and k, and r the 



2 



radial coordinate. Around typical neutrinospheres with p m io n ~ 12 g/cm 3 and the electron 
fraction per baryon Y e « 0.1, the order of magnitude of Eb is estimated to be 

Y e \—s ( p \ a / B \ 



which is indeed small so that a perturbation calculation can be applied. We caution the 
reader that the neutrino energy k as well as the electron number density N e depend on r. 
The deeper the neutrinos are emitted from, the higher their energies are. To the zeroth 
order, i.e. in the electron background without magnetic field, the resonance condition is 
satisfied at r$ as 

A 2 ,cos2# = 2V2G F N e (r )k(r ). (3) 
The presence of B field modifies the resonance position into 

r{4>) = r + S r cos0, (4) 

where 8 r <C r because Sb, the correction due to the presence of magnetic field, is very 
small(see Eq.(|D). Substituting Eq.@ into Eq.([l]) and taking into account the r-dependence, 
we obtain, to the leading order in eb or 5 r /r , 



dlniV; 



dr 



+ 



d In k 



dr 




Uj?' B - < 5 ' 



where we have used the fact that dN e /dr < and dk/dr < 0. Under the assumption 
of blackbody radiation for neutrinos, the energy flux of neutrinos T v {r) is proportional to 
T 4 so that |dln/c/dr| = 4|dlnT/dr|. Defining the scale heights at the resonance point as 
hr = (dlnT/dr]" 1 and h Ne = |d In N e /dr\~^ and using N e ps pI/Stt 2 where \i e is the chemical 
potential of the degenerate relativistic background electrons, Eq.(||) reduces to 



ZeB h^h 



(6) 



(7) 



2pi 2 e 4h Ne + h T 

The asymmetry of the total neutrino momentum distribution is, then, given by 

Ak 1 Jq T v (r + 5 r cos (p) cos (p sin i 
k 6 Jq T v (r + S r cos 0) sin (j. 

« 01 ( I ) ( B - ) (^-Y 

\l + (h T /4h Ne ) J V3 x 10 15 Gauss/ VllMeW ' 

where all six species of neutrinos are approximated to equally contribute to the total energy 
of a supernova. It is to be emphasized here that unity in the denominator of the first 
bracket in Eq. (^) is due to the r-dependence of neutrino energy. Note that we have not used 
the relation diV e /dT m (dN e /dT)^ e in Ref. [|IJ which is, strictly speaking, not valid inside 
protoneutron stars as pointed out in Ref. 0. The neglect of r-dependent neutrino energy, 
which amounts to replacing 1/(1 + (/iy/4/ijv e )) by 4hN e /hr, leads to 
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which is the result in Ref. M. Further, if we ignore the temperature dependence of chemical 
potential which justifies dN e /dT « (<9iV e /<9T) M( ., the ratio of the scale heights can be written 
as h Ne /h T ~ 3iV e /2T 2 /i e . Then, using the relation N e ~ /ig/37r 2 , we can rewrite Eq.(H) as 



Ak 
~k 



B 



2vr 2 T 2 



(9) 



p ~ 10 
Ak/k r 
are 



(with the correct flux factor 1/3). If we set B ~ 10 15 Gauss, 
1 g/cm >5 , Y e ~ 0.1(/i e ~ HMeV) and T ~ 3MeV, we have h^ e /hr ~ 0.7 and therefore 
^ 0.01. However, numerical supernova calculations |10| show that the scale heights 



which is the result in Ref. 

-ill,-. / 3 



6 km, h p ~ 7 km, and hx ~ 22 km at the stage of the protoneutron star in a 



of the scale heights of the model [10 



model of fixed baryon mass, 1.4M , without accretion or convection. Therefore, the ratio 

■> 0.3 which is slightly different from the 
3iV e /2T 2 /i e , increasing the value of B to 
. If hN e /hr turns out to be large 



is h N Jh T < 

result, 0.7, of the approximate relation h^Jhx ~ 
have Ak/k ~ 0.01. This was the point stressed in Ref 
(h^ e /hx > 1), then the r-dependence of neutrino energy becomes relatively important and 
one has to resort to a general formula in Eq.(^). In general, the value of h^Jh^ depends 
on the model of the birth of neutron star and on the time evolution of the collapse-driven 
supernova explosion. In summary, the incorporation of the r-dependence of the neutrino 
energy as manifested in the first parentheses in Eq.f^J) requires, in general, higher B values 
in order to explain the observed birth velocity of pulsars. 

Now, let us discuss the problem in a three-neutrino generation scheme with a mass 
hierarchy that can accommodate the solar neutrino data |llfl. According to a general theory 



of three neutrino oscillations, the weak eigenstates of neutrinos v e ,u, T are superpositions of 



the mass eigenstates ^1,2,3 via a 3 x 3 unitary matrix U [12]: 











= u 




V T ) 




V ^3 



(10) 



Assuming, for simplicity, that CP is conserved, we adopt a useful parameterization of U as 



U 



(11) 
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where = cosu, = sinu; and the A's are the Gell-Mann matrices. In an electron-rich 
medium, the equations of motion for the weak eigenstates of neutrinos are 



dt 



1 

2E 



U 



( ml 


V 







mt, 




(12) 



where A is the matter-induced mass squared of v e which is, in the presence of magnetic field, 




A = 2V2G F N e E + V2eG F — ^ k ■ B. (13) 



We will consider the pulsar kick velocity by three neutrino mixing under two natural 
conditions: (1) ni\ -C m2 -C m^, (2) the mixing angles in vacuum are small. It is easy 
to see then that three neutrino oscillations in matter separate into two parts, each being 
similar to two neutrino oscillations in matter. Two resonant conversions can take place at 
well separated locations. We call Rh the resonance surface which occurs at higher density, 
and Ri the one at lower density. The conditions for each resonance are [O] 



\ {m?(l + C u ) + m|(l - C W )}1 C 2v , (14) 
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ml - ml)C 2lJJ 



In a three generation scheme, therefore, the relative positions of two resonance surfaces 
to neutrinospheres are essential in order to study the effects of neutrino oscillations on 
the birth velocities of pulsars. First, let us discuss the relative positions of three neutri- 
nospheres. A neutrinosphere may be considered as an imaginary surface such that inside 
the surface neutrinos are trapped due to the interactions with medium, but outside they 
escape freely. Therefore, the larger the scattering cross section is, the outer the neutri- 
nosphere is located. Since the temperature (T pa 3 MeV) around the neutrinospheres is 
too low to thermally create \x or r, the neutrinosphere is expected to coincide with the 
v T neutrinosphere, S T . As can be seen in Eq. fll4D , the locations of resonance surfaces are 
dependent on the masses and the mixing angles of neutrinos in vacuum which have been 
restricted by various measurements and observations. According to a comprehensive analy- 
sis of solar and atmospheric neutrino experiments in a hierarchical three-generation scheme, 
mi "C TU2 <C TU3, the Mikheyev-Smirnov-Wolfenstein (MSW) solution to the solar neutrino 
deficit problem suggests m\ — m\ pa 10~ 5 eV 2 [fyj and the atmospheric neutrino anomaly 
leads to m\ > 10~ 3 eV 2 fT5 |. Therefore, the density around neutrinospheres with p pa io 11 ~ 12 



g/ cm 3 is too high for the lower resonance to occur, implying that the lower resonance surface 
Ri lies outside the S e . Within the current experimental limit of m 3 , the Rh (the resonant 
conversion v e «-> z/ T [J) can lie between the S e and the S T fl6fl . 



In the following we assume that the Ri for v e <-> is located outside the S e , but the 
Rh for v e *r* v T between the S e and the S T . At the Rh, v T produced inside the S T will be 
converted into u e and absorbed, whereas u e produced inside the S T will turn into v T and will 
propagate out with the collapsing energy of the protoneutron star. The does not undergo 
any resonant conversion at the Rh- As neutrinos pass the lower resonance surface, neutrino 
flavors are resonantly converted again on the magnetically distorted surface. However, this 
resonant conversion outside the S e does not contribute to the pulsar velocities because the 
momentum transfer of neutrinos is independent of whether or not the neutrino flavor is 
changed. This is in contrast to the case where the resonance surface lies between the S e 



Of course, this is an approximate interpretation, which is valid when the neutrino mixing angles 
in vacuum are very small. 
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and the S T so that v T converted from v e propagates out but v e from v T cannot, leading to 
an asymmetry in the neutrino momentum distribution. Once neutrinos escape from the S e , 
therefore, flavor-changing of the neutrinos does not affect the kick velocity of pulsars. We 
summarize that if the value of (m| — m 2 ) is restricted by the MSW solution to the solar 
neutrino problem so that the lower resonance occurs outside the S e , the consideration of 
three neutrino oscillations does not affect the result of the KS mechanism based on the two 
generation scheme. 

Next, we discuss the effects of possible non-adiabatic neutrino oscillations on the birth 
velocity of pulsars. It is, of course, obvious that non-adiabatic processes make u e —>■ v T 
conversion less efficient, reducing the kick velocity. We wish to study this in some detail. 
Let V{y a — > up) be a general probability for u a produced deep inside the resonance surface 
to be observed as up outside the resonance surface. In adiabatic cases where instantaneous 
effective mass eigenstates approximate the energy eigenstates of Hamiltonian, the transition 
probability VpSy a — > up) |]T7| is 



C 2 S 2 

V A (u a -> vp) = {vp\ { s { c { ] ( ^ \ ) ( ~* 2 ^ 2 ) \p a ), (15) 




where Cf^ = cosd f t i and 9f(9i) is the effective mixing angle of neutrinos in matter after 
(before) neutrinos go through the resonance surface and 



ke) = I J, k> = ( 1 ]■ (16) 

If the oscillation process is non-adiabatic, a level-crossing between the effective mass eigen- 



states takes place, particularly near the resonance point [12]. With the level-crossing 



(Landau-Zener) probability Plz-, we have the following non-adiabatic transition probabilities 

PnaK - u p ) = (vp\ ( ^ ^ ) ( 1 ~p P L L z Z ! P _ L p LZ ){ C S }c 2 ) l^>- ^ 

The level-crossing probability is Plz = exp(— jQ) where the adiabaticity parameter Q is 
given by Q = A 2 ] S% g h,N e /kC28 in the two generation case. For Aq ~ 10 4 eV 2 , hjsr e ~ 6 
km, and k ~ 10 MeV, the non-adiabatic effect cannot be ignored when the vacuum mixing 
angle 9 is smaller than 10~ 4 (Recall that the mixing angle of relevance is v e — v T mixing 
angle, therefore, a very small mixing angle is not ruled out). For example, the level-crossing 
probability is Plz ~ 0.38 for 9 = 10~ 4 . In the extremely non-adiabatic case, flavor changing 
processes would not take place, so that the KS mechanism becomes inoperative. 

Let us consider a general non-adiabatic case with a possibility that the fluxes of z/ e 's 
and z/ T 's from the protoneutron star are different. The collapsing supernova core has about 
10 57 protons which produce neutrinos via p + e~ — > n + u e , resulting in the formation of a 
neutron star. Since each v e emitted from the core carries away an average energy 10 MeV, 
about 10 52 ergs are emitted by i/ e 's during the neutronization processes. This is less than 
5% of all the u's radiated [|T^]. The rest of the neutrinos come from pair processes such as 
e + + e~ — > Ui + Ui, where i = e, \i or r. Since and v T are produced via neutral currents but 
v e via both charged and neutral currents, there exist some differences between the v e flux and 



6 



other neutrino fluxes, defined collectively by v@. We will introduce the parameter a Ve such 
that neutrino fluxes are given by T Ve = (1 + a Ue )J- and T Vj3 = J 7 , where T is the common 
flux factor. The existence of another enhancing channel for v e production implies a Ue > 0. 
The asymmetry in the neutrino momentum distribution due to non-adiabatic oscillations is 
then given by 

AAA ^ ^NApe -> Vt) (1 + CQ Jo F( r + $r COS 0) COS SU1 0d0 

k / NA (6 + a„J Jo 71 J^(r + S r cos 0) sin 0d0 




where the second approximation is valid for small value of a Ue . The term (Ak/k) on the right 
hand side of Eq.flTED is given by Eq.(|7|). In the small vacuum mixing case (Of ~ 0, 0{ ~ 7r/2) 
the transition probability T-Wpe — ► &v) can be further approximated, resulting in 




This shows that the non-adiabaticity has a tendency to decrease the asymmetry in the neu- 
trino momentum distribution. Equation (|19|) shows that extremely non-adiabatic processes 
(Plz = 1) cause no kick velocity due to the absence of flavor conversion, as discussed before. 
The pulsar velocity increases with a Ue . This is also expected because the more z/ e 's are 
produced inside the core, the more they turn to the ^ T 's, leading to a larger birth velocity of 
pulsar. Therefore, the effect of non-adiabaticity is opposite to that of a Ue > 0. When we set, 
for example, the vacuum mixing angle = 1.6 x 10 -4 and a — 10 %, the contributions of 
the non-adiabaticity and the excess flux of electron neutrinos to the momentum asymmetry 
are almost canceled out. This cancelation depends very sensitively on the value of 0. 

In summary, we have studied the Kusenko-Segre mechanism in a more general condition 
than in Refs. [0 and 0. First, the position dependence of neutrino energy has been taken 
into account in the calculation of the neutrino momentum asymmetry to explain the peculiar 
pulsar velocity. This aspect becomes relatively important in the case hN e /hr > 1- In 
addition, we have clarified the difference between the two treatments in Refs. Q and |§. 
The KS mechanism is also studied in a three generation scheme and in a non-adiabatic case. 
If the solar neutrino and atmospheric neutrino problems are to be solved with the standard 
mass hierarchy, a generalization of the treatment in Ref. in a three-generation scheme 
leads to no new modifications. The non-adiabaticity has a general tendency to decrease the 
kick velocity whereas a larger electron neutrino flux compared to that of other flavor gives 
rise to a higher momentum asymmetry. 
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